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Abstract - Recurrence networks are a novel tool of nonlinear time series analysis allowing the 
characterisation of higher-order geometric properties of complex dynamical systems based on 
recurrences in phase space, which are a fundamental concept in classical mechanics. In this Letter, 
we demonstrate that recurrence networks obtained from various deterministic model systems as 
well as experimental data naturally display power-law degree distributions with scaling exponents 
7 that can be derived exclusively from the systems' invariant densities. For one-dimensional maps, 
we show analytically that 7 is not related to the fractal dimension. For continuous systems, we 
find two distinct types of behaviour: power-laws with an exponent 7 depending on a suitable 
notion of local dimension, and such with fixed 7 = 1. 



Introduction. — Power-law distributions have been 
widely observed in diverse fields such as seismology, econ- 
omy, and finance in the context of critical phenom- 
ena [l]-[3]. In many cases, the underlying complex systems 
can be regarded as networks of mutually interacting sub- 
systems with a complex structural organisation. Specifi- 
cally, numerous examples have been found for hierarchical 
structures in the connectivity of such complex networks, 
i.e., the presence of scale-free distributions P{k) ^ k~^ of 
the node degrees [4||5]. Such hierarchical organisation is 
particularly well expressed in network of networks, or in- 
terdependent networks, which constitute an emerging and 
important new field of complex network research [6|u]. 
The interrelationships between the non-trivial structural 
properties of complex networks and the resulting dynam- 
ics of the mutually interacting subsystems are subject of 
intensive research [8J[9]. 

Among other developments, one of the main recent 
achievements of complex network theory are various con- 
ceptionally different approaches for statistically charac- 
terising dynamical systems by graph-theoretical meth- 



ods 10 ■ 13 . In this Letter, we report and thoroughly 



explain the emergence of power-laws in the degree dis- 



tribution of so-called recurrence networks (RNs) [13 16 
for various paradigmatic model systems as well as exper- 
imental data. RNs encode the underlying system's re- 
currences in phase space and are based on a fundamen- 
tal concept in classical physics 17 . Due to their direct 



link to dynamical systems theory, RNs are probably the 
most widely applicable type of complex networks inferred 
from time series introduced so far. Although the system's 
temporal evolution cannot be reconstructed from the RN, 
this representation allows for an analysis of the attractor's 
geometry in phase space using techniques from network 
theory. Specifically, nodes represent individual state vec- 
tors, and pairs of nodes are linked when they are mutually 
closer than some threshold distance e > flS] (cf. Fig.fTl). 
According to this definition, RNs are random geometric 
graphs '19] (i.e., undirected spatial networks [20| ), where 
the spatial distribution of nodes is completely determined 
by the probability density function of the invariant mea- 
sure of the dynamical system under study, and links are 
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Fig. 1: (Colour online) Construction of a RN: (A) time series 
xt, (B) reconstructed two-dimensional phase space trajectory 
with state vectors yt = {xt,Xt+i) and an exemplary e-ball 
around y2 (dashed) for e = 3. (C) Resulting recurrence net- 
work. 



established according to the distance in phase space. Con- 
sequently, their degree distribution P{k) directly relates to 
the system's invariant density p{x). 

In this work, we demonstrate the emergence of scaling in 
the degree distributions of RNs and provide some evidence 
that this phenomenon is (unlike many other scaling expo- 
nents occurring in the context of dynamical systems) com- 
monly unrelated to the fractal attractor dimension, except 
for some interesting special cases. Instead, the power-laws 
naturally arise from the variability of the invariant density 
p{x) of the system {i.e., peaks or singularities oi p), as we 
will show numerically as well as explain theoretically. 

Power-law scaling and singularities of the invari- 
ant density. — As initial examples, Fig. [2] illustrates 
the presence of power-law degree distributions in the RNs 
obtained for several prototypical low-dimensional chaotic 
systems with a suitable choice of the systems' character- 
istic parameters: (i) the Rossler system in spiral-chaos 
regime: i — — y— z, y = x + 0.2y, z — {i.2 + z(x — 5.7); (ii) 
the Lorenz system: x = 10{y — x), y = x{r — z) ~ y, z ~ 
xy~8/3z; and (iii) the Henon map: Xn+i = 1 — l-4a;^-|-y„, 
Vn+i = 0.3a;„. For the Rossler and Lorenz systems, a 
proper time discretisation has been used as explained in 
detail in the figure caption. 

In all cases, scaling emerges only if the distance thresh- 
old e is chosen small enough, which corresponds to a small 
average degree (k) and hnk density p = (k)/{N— 1) of the 
resulting network (N being the number of nodes) . We note 
that the respective range of e should be sufficiently higher 
than the threshold for which a giant component exists [25] . 
The size of the scaling regime decreases with growing p 
and becomes hardly detectable for p > 1% (Fig. [2]). The 
distance threshold e also occurs in dimension estimation 
where the limit e — > is taken {e.g., 26 ). In contrast, a 
RN is based on one finite e. Our smallest e are, however, 
still large enough to avoid the problem of lack of neigh- 
bours 



23 since the correlation integral C{e) still shows 



10" 



10"' 
10"' 




'% 


C» + 


/D,.1.92 


10" 


10-= 



Bio° 




10 10 10 10 10 10 10 10 

k k 

Dio° 



O lU 






-»»,— 




£io-^ 


10-' 
10-^ 


^^ Dj.2.03 




in-* 




10^ 


10-' 1C 


~ 



^o'■ 



10" 



10' 



10' 



10" 





-^. "X 














10-' 




X 


%Vr. + 


10-' 


x' 


D^ -1 .20 






10-* 


10-' 1C 


-' 



10"' 



10" 



10' 



10' 



Fig. 2: Cumulative degree distributions -F(fc) — J^^^j. -P(fe') of 
the RNs obtained for several discrete maps: (A) a;-component 
of the first return map of the Rossler system (7 = 2.16 ± 0.03), 
(B) map of consecutive local maxima of the z-component of 
the Lorenz system with r = 28 (no scaling), (C) as in (B) with 
r = 90 (7 = 2.64±0.18), and (D) Henon map (7 = 2.88±0.04). 
Estimates of 7 > 1 have been obtained from the cumulative 
degree distributions F(k) ~ k} ' (in the following figures, we 
show either -F(fc) or P(fe) when sufficiently straight in a log- 
log plot) by means of a maximum likelihood approach pi] as 
averages considering 100 different values of p for which a power- 
law appears. For each p, 5 different realisations with random 
initial conditions are used. We consider RNs of size N = 2x 10^ 
and the Euclidean norm in all examples discussed in this Letter. 
Four cases were chosen for illustration corresponding to a link 
density of pi = 0.02% (o), pa = 0.03% (•), pa = 0.05% (>), 
and pi — 3% (+). A power-law is hardly detectable for p4. 
Insets: corresponding log-log plot of the correlation sum C(e) 
vs. e, where the correlation dimension D2 is estimated by linear 
regression and agrees well with values from the literature |22| - 



in Fig.[2|. 

A theoretical explanation of the emergence of power-law 
distributions of RN-based degree is based on the general 
theory of random geometric graphs 
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where nodes are 
sampled from some probability density function p{x). For 
a RN, the space is the phase space of a dynamical sys- 
tem and the nodes are states sampled at discrete times. 
If we assume that the system is ergodic, the sampled tra- 
jectory is already close to its attractor, and the sampling 
times are generic (particularly, the sampling interval is 
co-prime to any period lengths of the system), the nodes 
can be interpreted as being sampled from the probability 
density function p{x) of the invariant measure p, of the at- 
tractor 
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The degree distribution of a general random 
geometric graph, P{k), is derived from p{x) in the limit of 
large network size N as fl9| 



the same dependency on e as for larger values (see insets 



P{k) = / dxp{x)e-°'P^''\ap{x)f/k\ 



(1) 
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with a = (k) / J dx p{x)'^ (note that the computation of 
(fc) involves integration of p{x) over the e-neighbourhood 
of all points x and thus implicitly depends on the specif- 
ically chosen e as well as the sample size N). Hence, the 
invariant density of the system exclusively determines the 
existence of a power-law in P{k) and its exponent 7. 

For systems with a one- dimensional phase space, it can 
be shown that under some weak conditions on p{x) it 
holds: 

Pik)^'^E.er>-H'^)\p'i^r'- (2) 

This implies that if p{x) has a power-law-shaped peak at 
some state xq, i-e., p{x) ~ jx — xo\~^''^ for some 7 > 0, 
the degree distribution P{k) also follows a power-law but 
with the reciprocal exponent, P{k) ~ k~'^ . Specifically, a 
slower decaying invariant density leads to a faster decaying 
degree distribution. Note that not all invariant densities 
lead to power-laws: If p{x) is Gaussian, we get P{k) « 
2(— 21n/c-\/27r/a)~^/^ instead of a power-law. 

More generally, we can deduce that the presence of sin- 
gularities in the invariant density is the key feature deter- 
mining whether or not the resulting RN has a power-law 
degree distribution. This relationship can be intuitively 
understood: If p{x) has a singularity at some point xo in 
phase space, then a time series of the associated dynam- 
ical system will return very often to the neighbourhood 
of Xq. Hence, nodes with a high degree will accumulate 
close to the singularity. If the resulting invariant density 
obeys a power-law decay, Eq. ^ implies the emergence of 
a power-law degree distribution. If there is more than one 
singularity oip(x), one can expect the resulting degree dis- 
tribution being related to a weighted sum of the influences 
of these points. Vice versa, the presence of a power-law de- 
gree distribution in the RN of a dynamical system requires 
the existence of a power-law in the invariant density, i.e., 
the presence of a singularity. We therefore conjecture that 
a local power-law in the invariant density is a necessary 
and sufficient condition for the emergence of a scale-free 
RN. Beyond the explicit results for one-dimensional sys- 
tems as discussed above, we further show below that the 
emergence of a power-law scaling is also possible in higher 
dimensions and conjecture that this requires the presence 
of a dynamically invariant object {e.g., an unstable or hy- 
perbolic fixed point) close to which the invariant density 
scales as a power-law at least in one direction. 

Power-law scaling vs. fractal dimension. As an 

example for discrete-time systems leading to power-laws, 
consider the generalised logistic map f28J 



f{x) = l^\2x~\ 



P 



(3) 



with /3 ^ |. For /3 = i, 1, 2, this gives the cusp map, tent 
map, and standard logistic map, respectively. For general 
/3 > 0, the unit interval [0, 1] is mapped onto itself by a 
symmetric function with a maximum of 1 at a; = ^, thus 
having two pre-images for each a:; < 1. For /3 > 1, the 
associated invariant density p{x) has two peaks &t x — Q 
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Fig. 3: (Colour online) (A): Cumulative degree distribution 
F{k) of the RNs for the map ^ with /3 = 2, 2.5, 3, 4. Note that 
we divide the degree fc by e so that the a;-axis is e-dimensionless. 
(B) P{k) of the cusp map P — 0.5. The inset shows the corre- 
sponding iterative function f{x). The dashed line of (B) has a 
slope of 1. The legend indicates the respective link densities p. 



and X = 1 with p{5) = p{l - S) ^ ^(i-/^)//? for small 5. 
Hence, the degree distribution P{k) shows a power-law 
with the exponent 



7 = /3/(/3-l). 



(4) 



Numerical results shown in Fig.|3]4. for several different 
values of (3 agree precisely with Eq. Q. In contrast, 
for /3 = 1 the nodes are uniformly distributed, and the 
degree distributio n d erived from Eq. (fTl) is Poissonian, 



P(fc) 



-"a'=/A:! 
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For /3 = 5 , we get p{x) — 2 ~ 2x 
[SOJfSl], which lea(£s~to a specific type of "power-law" in 
Piji) with 7 = — 1 as shown in Fig.lsb. These results im- 
ply that the scaling exponent is not simply related to the 
fractal dimension: the attractor has the box-counting di- 
mension Dq = 1 independently of /3, whereas 7 changes 
with varying /3 (Eq. Q). However, the correlation di- 
mension Di also depends on /3 (Z?2 = 1 for ^ < 2, and 
D2 = 2//3 if /3 > 2 [23]), i.e., there is an indirect rela- 
tionship between 7 and D-i for certain special cases. The 
different behaviour of the mentioned dimensions results 
from the fact that Dq exclusively considers the number 
of boxes required for covering the attractor, but not their 
individual probability masses as Di and other notions of 
fractal dimensions do. 

Turning to continuous-time systems, we next compare 
the above findings with those for some dis cretis ed standaid 
examples. On the one hand, for the Rossler system, we 
consider the successive a;-values when passing the Poincare 
section at y — with y < 0. As shown in the inset of 
Fig.[4K, the resulting first return map has a shape similar 
to the case of /3 = 1.87 in Eq. ([3|. Hence, we expect a 
power-law with the exponent 7 sa 1.87/(1.87 — 1) = 2.15. 
The invariant density has several dominant peaks, which 
are together responsible for the power-law observed in 
Fig.[2K with 7 indeed close to 2.15. In fact, P{k) is a 
mixture of individual power-laws corresponding to the in- 
dividual peaks of p{x), whose exponents are all roughly 
the same. On the other hand, for the Lorenz system, 
we obtain a one-dimensional map by studying the local 
maxima z^^^ of z for successive cycles 32 , i.e., mapping 
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Fig. 5; (Colour online) (A): P{k) for the RN of the Rossler sys- 
tem with screw-type chaos 34*. Three different link densities 
are chosen for illustration as indicated by the legend, a larger 
link density p4 does not produce a power-law (inset). The 
slope of the dashed line is —1.33. (B): t?-capacity dimension 
Do ~ 1-30 in three small cuboidal neighbourhoods of differ- 
ent size (in terms of phase space distance in each coordinate 
direction, see the legend). 



Fig. 4: (Colour online) PDF of (A) the successive x values 
passing the plane oi y = of the Rossler attractor, (B,C) PDF 
of Zmax for the Lorenz system of r = 28 (B) and r — 90 (C), 
and (D) marginal density of the x component of the Henon 
map. The insets show the associated first return plot. The 
dashed curve in (B) shows a fit by a Gaussian distribution. 



zlax to zl^l (inset of Fig.|4^). For r = 28, this first re- 
turn map has a similar shape as Eq. ([s]) for /3 = 0.5 (inset 
of Fig. [3J3), but the corresponding density is beh-shaped 
without a peak. Indeed, we do not observe a power-law for 
P{k) in this case (Fig.[2j3). However, increasing r changes 
the shape of p{x) qualitatively. For example, at r = 90 
(Fig.|4]C) the density has peaks at several points, explain- 
ing the observed power-law in Fig.[2p. A similar behaviour 
can be observed for the Henon map, though the marginal 
invariant density of the a;-component has a more complex 
structure (Fig.jip). 

While there is no unambiguous relationship between 
7 and the fractal dimension already for discrete sys- 
tems, the situation becomes even more complicated for 
continuous-time systems which are not discretised via a 
Poincare section or otherwise. For two-dimensional flows 
X = $(a;) with only one peak in p{x), the respective type 
of behaviour depends on the eigenvalues of the Jacobian 
-D'l>(a;o) at the fixed point Xq as well as on the shape of 
p{x). Specifically, in many cases (that shall not be further 
discussed here) the existence of a power-law for P{k) can- 
not be evaluated easily, whereas in other cases, one can 
analytically derive a power-law with a very small expo- 
nent of 7 = 1. In turn, the following numerical results 
suggest that there are also examples displaying a distinct 
relationship between 7 and some suitably defined local di- 
mension: 

For the Rossler system in the regime of screw-type 
chaos with a homoclinic point at the origin fulfilling the 
Shilnikov condition 



33 34 



the invariant density is dom- 
inated by its peak at the origin. The degree distribution 
P{k) of the corresponding RN shows a power-law with 



7 w 1.33, which agrees fairly well with the ^-capacity 

(Fig. [5]). We also observe 



dimension Dq 



defined in 
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similar scaling laws for both numerical model and exper- 
imental data (output intensities) of a single-mode CO2 
laser 
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The underlying system has a saddle-focus S 
embedded in the chaotic attractor (Fig. |6]4.) which causes 
a spiking dynamics |36lf37]. The attractor is dominated 
by a homoclinic orbit emerging from and converging to S. 
The degree distributions P{k) resulting from both model 
and experimental data suggest power-laws with 7 sa 1.35 
(Fig. lob), which qualitatively agrees well with the point- 
wise dimension of the attractor around S. Finally, simi- 
lar results can be obtained for a predator-prey food-chain 
model with four competing species 
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which also dis- 
plays homoclinic chaos, where we observe 7 « 1.9 in agree- 
ment with Dq (Fig. w\. This variety of examples under- 
lines the general importance and wide applicability of our 
findings. 

Technical aspects. — In general, we have to make 
two cautionary notes on the numerical study of scaling 
laws in RNs. 

First, in the previous continuous-time examples, the 
presence of power-laws with the numerically estimated ex- 
ponents (see above) cannot be rejected on a 90% signif- 
icance level using Kolmogorov-Smirnov tests. However, 
the alternative of a power-law with 7=1 can also not be 
rejected at the same level. Hence, power-laws with quali- 
tatively different exponents describe the data comparably 
well. It remains an open problem to determine the correct 
7. We note that this is a general problem when evaluating 



hypothetical power-laws from finite data 21 39 



Second, experimental data often consist of only one 
measured variable. Hence, a reconstruction of the associ- 
ated phase space trajectory is necessary prior to RN anal- 
ysis, e.g., by time-delay embedding 40 . Like estimates 



of dynamical invariants or complexity measures 41 42 



the power-law behaviour of P{k) can depend on the par- 
ticular observable, because different coordinates of a dy- 
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Fig. 6: (Colour online) Experimental laser data: (A) Phase 
portrait in subspace of (xi, 2:2, aje), where the saddle focus 5* is 
at the most dense region. (B) P{k) of RNs from model data 
for three different link densities (see the legend); a larger link 
density p4 does not lead to a power-law (inset). (C) P{k) of 
RNs from experimental data (inset; cumulative distribution 
F{k)). (D) Point-wise dimension Df — 1.35. All dashed lines 
have slope —1.35. 



namical system often have different marginal densities. 
Specifically, embedding theorems ensure topological invari- 
ance {i.e., properties of the dynamical system that do not 
change under smooth coordinate transformations are pre- 
served), but no metric invariance of the attractor's geom- 
etry including p{x). For example, the logistic map and the 
tent map (/3 = 2, 1 in Eq. ([3|, respectively) are topologi- 
cally equivalent under the transformation x i-> sin {7tx/2), 
but the different invariant densities with respect to their 
original coordinates (that have been used for constructing 
the RNs from metric distances in their respective phase 
spaces) lead to distinct scaling exponents 7 (Fig.lsl 



Conclusions. — In summary, we have reported an in- 
teresting novel aspect of the geometrical organisation un- 
derlying the dynamics of many complex systems in physics 
and beyond. Specifically, we have provided an analytical 
explanation of the emergence of power-laws in recurrence 
networks constructed from sampled time series based on 
the theory of random geometric graphs. Unlike for compa- 
rable complex network approaches [lOJ[TI], this scaling is 
not simply related to the system's fractal dimension, but 
determined by both the singularities of the invariant den- 
sity and the considered spatial scale e. We emphasise that 
dimensions are defined in the limit of e — >■ and practically 
estimated by a series of e values, whereas the power-law 
exponent 7 of the RN appears for each sufficiently small 
e individually. Note that in contrast to the degree, the 
transitivity properties of RNs have a direct relationship 
with attractor dimension 
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In comparison with the invariant density itself, fractal 
dimensions are a rather specific characteristic. In particu- 
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Fig. 7: (Colour online) As in Fig. [5] for the Lotka-Volterra 
system [38] with homoclinic chaos. The slope of the dashed 
line in panel (A) is —1.90. 



lar, they do not simply describe the whole system (as the 
invariant density itself does), but quantify density vari- 
ations on the attractor viewed at different spatial scales 
e [26II44]. Conversely, the scaling exponent 7 directly 
characterises a power-law decay of the density in phase 
space independent of a specific scale. In this spirit, both 
fractal dimension and scaling exponent 7 capture concep- 
tually different aspects of the geometric organisation of a 
dynamical system in its phase space. However, although 
there is no general relationship between 7 and fractality, in 
some special cases the power-law exponent coincides with 
some notion of dimension. This has been demonstrated 
for several example systems as well as experimental data. 
In turn, we have found that in other cases the value of 7 
drops to 1. Further studies are necessary in order to better 
understand this complex relationship between power-law 
degree distributions and fractal scaling (i.e. under which 
general conditions related to a system's structural organ- 
isation both scaling exponent and fractal dimension coin- 
cide), particularly in continuous dynamical systems. 

From a conceptual perspective, we would like to remark 
that studying a single scalar property like the scaling ex- 
ponent of a recurrence network or the fractal dimension 
cannot provide a complete view on the structural organ- 
isation of a nonlinear complex system. Specifically, both 
characteristics capture distinct and complementary fea- 
tures related to the probability density of the invariant 
measure. In this spirit, the power-law exponent 7 quan- 
tifies a fundamental property that has not been explicitly 
studied so far. Because its relationship with the features 
of possible singularities of the invariant density is intuitive 
(i.e. the emergence of power-law degree distributions has 
some clear physical meaning), one particular strength of 
studying the degree distribution of recurrence networks is 
that it potentially allows identifying the presence of such 
singularities in complex situations {e.g., for observational 
data). 
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